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ABSTRACT 

The  Air  Force  Research  Laboratory  Munitions  Directorate  (AFRL/RW)  is  performing  a  series  of  studies  to  validate 
the  ability  of  codes  to  accurately  predict  the  impulse  response  of  simple  mechanical  structures,  i.e.,  impact-loaded 
bars  and  plates,  over  a  wide  frequency  range.  The  modal  responses  of  both  structures  are  predicted  using  finite 
element  codes  and  simple  analytic  calculations.  An  experiment  on  isolated  (i.e.,  free-free)  structures  using 
impulse  hammers  as  the  excitation  is  described.  The  outputs  of  high  bandwidth  strain  gages  and  accelerometers 
monitor  the  elastic  wave  propagation  in  the  structures.  The  corresponding  frequency-domain  characteristics  (e.g., 
the  power  spectral  density  of  each  sensor  as  well  as  the  transfer  functions  between  them)  are  measured  over  a 
wide  frequency  range  in  a  cylindrical  bar.  Leakage  is  found  to  significantly  reduce  sensor  coherence  over  a  broad 
spectral  band.  Use  of  an  exponential  window  restores  the  coherence  in  the  original  signal,  leading  to  a 
broadband  coherence  below  10  kHz,  but  this  approach  broadens  the  modal  peaks  observed  and  introduces 
errors  in  damping  estimates.  The  initial  results  from  an  ongoing  effort  to  characterize  the  modal  response  of  a 
hanging  plate  are  also  presented. 

INTRODUCTION 

Although  modeling  and  simulation  (M&S)  tools  are  pervasive  for  mechanical  design,  the  ability  of  M&S  tools  to 
accurately  predict  impact  loads  remains  questionable.  For  example,  explicit  finite  element  shock  physics  codes 
are  not  well-characterized  for  high-frequency  coherence.  Recent  studies  [1,  2]  have  focused  on  quantifying  the 
long-duration  (i.e.,  many-wavelength)  propagation  of  impulses  in  simulations;  initial  results  have  shown  that 
underlying  numerical  filtering  and  roundoff  errors  cause  waves  to  lose  coherence  rapidly,  especially  at  high 
frequencies.  This  represents  a  fundamental  limitation  to  predictive  analysis  of  impulsive  loading  and/or  the 
prediction  of  high-frequency  local  resonances  in  structural  assemblies. 

This  effort  represents  the  initial  results  from  ongoing  experimental  studies  used  to  validate  the  efforts  of  Idesman 
et  al.  [3].  Subsequently,  the  objective  is  to  provide  measurements  of  the  vibration  and  transient  response  in 
impact-loaded  structures,  such  as  simple  bars  and  plates,  from  impact  loads.  The  elastic  stress  wave 
propagation  is  characterized  using  an  impulse  hammer,  strain  gages,  and  accelerometers  with  high  bandwidth 
measurement  techniques  (fc  »  1 0  kHz);  the  output  is  characterized  with  respect  to  modal  (frequency)  as  well  as 
transient  responses  (e.g.,  the  dispersion  of  the  waves). 

THEORY 

Vibrational  Modes 

Two  systems  were  studied,  a  Ti  alloy  bar  (15.2  cm  long,  2.54  cm  diameter)  and  an  Al  alloy  plate  (61  cm  x  122  cm 
area,  0.6  cm  thick).  The  bar  is  ASTM  Grade  5  (6%  Al,  4%  V)  Ti  alloy  ( p  =  4.43  g/cm3,  £  =  114  GPa,  v=  0.33)  and 
the  plate  is  2024  Al  alloy  (p  =  2.73  g/cm3,  E  =  72  GPa,  v  =  0.34). 


For  a  simple  free-free  beam  element,  the  bending  modes  can  be  calculated  from  Harris  [4]; 
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where  an  is  a  constant  for  the  r?th  mode  (where  an  =  22.4,  61.7,  and  121  for  the  first  three  modes);  E  is  the  elastic 
modulus;  /  =  7zd4 / 64  is  the  area  moment  of  inertia  for  a  right  circular  cylinder  of  diameter  d;  ju  =  ^4/?  is  the 

mass  per  unit  length,  and  L  is  the  length  of  the  bar.  The  corresponding  predictions  are  4.97  kHz,  16.7  kHz,  and 
26.9  kHz  for  the  first  three  bending  modes. 

The  first  three  bending  modes  calculated  by  an  LS-DYNA  eigenvalue  analysis  (using  a  coarse  mesh  with  3350 
nodes  and  2646  elements)  are  4.57  kHz,  11.4  kHz,  and  20.0  kHz,  shown  in  Figure  1  below.  The  first  three 
extensional  modes  are  16.6  kHz,  33.0  kHz,  and  48.9  kHz.  Simulations  of  the  free-free  plate  were  also  calculated 
via  LS-DYNA  eigenvalue  analysis;  the  first  ten  mode  shapes  are  presented  in  Figure  2  below. 

Free-free  bar  modes 


Figure  1.  Model  predictions  of  the  modal  response  of  the  free-free  bar. 


The  longitudinal  vibration  modes  of  a  free-free  bar  can  be  readily  calculated  using  a  one-dimensional 
approximation  [5]  from  the  longitudinal  elastic  wave  speed  (c)  in  a  bar, 

c  =  Je7^>,  (2) 

where  p  is  the  density  of  the  bar  material.  Equation  (2)  gives  c  =  5072  m/s  for  the  Ti  alloy;  this  gives  a  traversal 
time  of  At=30  ps  for  a  15.2  cm  bar.  The  corresponding  fundamental  frequency  is 


which  gives  16.7  kHz  and  33.4  kHz  for  the  first  two  longitudinal  modes,  closely  matching  the  simulation  results. 

Free-free  plate  modes 
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Figure  2.  Model  predictions  of  the  modal  response  of  the  free-free  plate. 


Transient  Analysis 

The  transfer  function  is  often  used  to  characterize  the  propagation  of  elastic  waves.  The  transfer  function  (also 
called  the  frequency  response  function  or  FRF)  is  the  frequency-domain  relationship  between  the  input  and 
outputs  of  a  specified  system  and  is  defined  as 


H(cd)  = 
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where  G(co)  is  the  transfer  function  and  X(co)  and  Y(co)  are  the  frequency-domain  input  and  output,  respectively. 
The  transfer  function  provides  information  on  the  relative  amplitude  of  the  signals  in  frequency.  Coherence  is  the 
ratio  of  the  cross  power  spectral  density  (Pxy)  with  the  auto  spectral  densities  (Pxx  and  Pyy)  of  the  input  and  output 
signals  [6],  i.e., 
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The  coherence  is  generally  used  as  an  indication  of  the  quality  of  the  data  [7],  and  can  be  interpreted  as  a  metric 
to  determine  if  signals  are  indeed  correlated  (i.e.,  all  energy  is  propagated  from  x  to  y)  without  spurious  energy 
sources  or  sinks. 


BAR  IMPACT  EXPERIMENT,  RESULTS,  AND  DISCUSSION 

The  Hopkinson  split  bar  technique  is  well  established  and  widely  utilized  for  research  of  shock  waves  [8],  The 
experimental  apparatus  is  a  shortened  Hopkinson  pressure  bar;  the  aspect  ratio  is  important  in  this  validation 
experiment  because  it  yields  finite  element  simulations  with  fewer  elements.  The  excitation  for  the  bar  is  provided 
via  an  ICP  impulse  hammer  (PCB  Model  086C01,  calibration  factor  12.43  mV/N)  with  a  metal  tip  [9]  at  the  center 
of  the  circular  face  end.  To  enhance  precision  and  accuracy  of  the  impact,  the  impulse  hammer  was  suspended  in 
plastic  rings  to  create  an  impact  pendulum.  The  hammer  is  raised  up  to  30  degrees  (chosen  to  avoid  channel 
saturation)  and  released.  Figure  3  depicts  the  support  system  of  both  the  impulse  hammer  and  cylindrical  bar; 
the  plastic  rings  and  rubber  bands  minimizes  support  interaction  and  approximates  a  free-free  boundary 
condition. 
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Figure  3.  (a.)  Schematic  of  the  experimental  setup  and  support  system,  (b.)  Picture  of  the  “pendulum”  impact 

hammer  with  (c.)  detail  showing  center  alignment. 


The  bar  was  highly  instrumented  with  sensors  to  capture  the  elastic  wave  propagation  and  subsequent  vibration 
response.  Two  types  of  sensors  were  used:  semiconductor  strain  gages  and  piezoresistive  accelerometers.  An 
Endevco  model  7270A-2k  [10]  piezoresistive  shock  accelerometer  (2  kg  full  scale,  6.971  pV/m/s2  sensitivity, 
90  kHz  bandwidth)  was  bolted  on  the  free  (non-impacted)  end  of  the  bar.  The  side  of  the  cylindrical  bar  was 
instrumented  with  five  semiconductor  strain  gages,  positioned  1  in.  apart  from  each  other  and  the  bar  ends;  the 
strain  gages  and  their  labels  are  shown  in  Figure  4(a.).  Epoxy  was  used  to  mount  the  strain  gages  onto  the  rod 
(see  Figure  4(b.))  and  then  allowed  to  cure  for  24  hours.  The  leads  were  soldered  to  enamel  wire  and  electrical 
tape  used  to  ensure  there  was  no  electrical  contact  between  the  wires  and  test  object.  The  particular 
semiconductor  strain  gages  used  are  Kulite  AFP-500-090  bare  silicon  [11]  which  have  a  gage  factor  of  140  ±  5% 
at  75  °F  and  have  a  nominal  resistance  of  500  Q.  The  rise  time  of  the  semiconductor  sensors  is  approximately 
100  ns,  yielding  an  effective  bandwidth  of  5  MHz. 


(a.)  (b.) 


Figure  4.  (a.)  Picture  of  the  instrumented  bar  under  test  (the  impact  face  is  on  the  bottom  right)  and  (b.) 
schematic  of  the  attachment  scheme  for  the  semiconductor  strain  gages. 


Typical  data  from  the  impact  testing  is  shown  in  Figure  5  below;  the  data  is  shown  for  one  test  from  the 
accelerometer,  strain  gages,  and  the  force  hammer.  The  data  is  filtered  at  50  kHz  to  aid  visualization.  The  time 
domain  traces  show  a  distinct  time-of-arrival  delay  relative  to  the  initial  elastic  wave.  The  inset  of  Figure  5  shows 
that  a  quick  estimate  of  the  delay  is  approximately  32  ps,  which  is  in  close  agreement  with  the  calculation  from 
Eq.  (2).  The  modal  response  of  the  structure  is  also  immediately  obvious.  Another  obvious  feature  from  the  time 
history  is  that  strain  gage  3  has  the  largest  amplitude  due  its  center  placement  on  the  bar;  it  experiences 
maximum  bending  for  the  first  two  modes  whereas  strain  gages  1  and  5  have  much  lower  relative  displacement. 


The  frequency  response  of  the  structure  is  shown  as  a  power  spectral  density  of  the  sensor  output  in  Figure  6 
(using  a  Matlab  routine  based  on  a  spectral  estimating  technique  by  Welch  [12]).  The  mode  structure  is 
immediately  obvious  with  the  first  several  modes  easily  discernible.  These  values  match  closely  to  the  analytic 
and  computational  mode  predictions  (see  Figure  1).  It  is  also  important  to  note  the  bandwidth  of  the  sensors  is 
only  due  to  the  global  roll-off  from  the  low-pass  filtering  at  a  corner  frequency  of  50  kHz.  The  only  exception  is  the 
impulse  hammer,  which  has  a  bandwidth  of  approximately  9.5  kHz  with  a  metal  tip.  However,  leakage  is  apparent 
in  the  broadened  peaks. 
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Figure  5.  Sensor  output  (filtered  at  50  kHz)  from  the  force  hammer,  strain  gages,  and  accelerometer.  The  time- 
of-arrival  of  the  elastic  wave  (inset)  is  approximately  32  //s,  which  agrees  with  the  predicted  time  of  30  //s. 


Power  Spectral  Density 


Figure  6.  Estimated  power  spectral  density  of  the  various  channels  from  the  entire  time  history.  All  data  is 

filtered  at  50  kHz. 


The  transfer  functions  were  initially  determined  between  strain  gages  5  to  3,  3  to  1,  and  5  to  1  to  evaluate  the 
quality  of  the  experimental  data;  the  data  is  shown  in  Figure  7.  The  transfer  functions  exhibited  a  strong  relative 
amplification  or  attenuation  near  the  principle  modes  as  would  be  expected  from  the  particular  bending  and 
extensional  mode  shapes  and  the  sensors  different  locations  on  the  bar.  However,  examining  the  coherence 
between  the  measurements  in  Figure  8  show  that  the  coherence  is  limited  to  extremely  narrow  frequency  bands 
with  coherence  ~1,  supporting  the  earlier  observation  that  leakage  is  dominating  the  signals. 
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Figure  7.  Transfer  function  estimate  between  the  strain  gages  (data  is  filtered  at  50  kHz). 
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Figure  8.  Mean -squared  coherence  of  the  data  in  Figure  5  (data  is  filtered  at  50  kHz). 


Windowing  is  typically  performed  on  sampled  data  to  reduce  the  effect  of  leakage  when  the  FFT  is  used  on 
discrete  data  where  the  time  history  does  not  attenuate  to  a  negligible  in  the  total  sample  time.  The  time  history 
from  this  weakly-damped  system  has  these  characteristics  and  is  windowed  using  an  exponential  (Gaussian) 
decay  window  [13]  can  be  applied  to  the  response  to  attenuate  the  signal  at  the  signal  beginning  and  end  in  the 
time  domain.  This  introduces  a  Gaussian  convolution  in  the  frequency  domain  and  correspondingly  broadens  the 
spectral  response  and  makes  it  difficult  to  estimate  of  the  dampening  for  the  modes.  The  results  of  the  windowing 
(using  an  exponential  decay  time  constant  of  450  /js)  on  the  coherence  are  shown  in  Figure  9  below.  It  is  evident 
that  all  of  the  measurements — impulse  hammer,  strain  gage,  and  accelerometer — are  coherent  below  10  kHz. 
The  filtering  and  higher  modal  density  makes  the  coherence  difficult  to  evaluate  above  10  kHz;  this  is  the  subject 
of  continuing  studies  and  is  beyond  the  scope  of  this  paper. 
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Figure  9.  Mean-squared  coherence  of  the  data  in  Figure  5  (filtered  at  50  kHz)  using  an  exponential  window. 


PLATE  IMPACT  RESULTS,  ANALYSIS,  AND  DISCUSSION 

In  addition  to  the  short  bar  impact  tests,  similar  analysis  was  performed  on  a  hanging  plate.  The  plate  was 
suspended  using  elastic  cables  and  instrumented  with  several  Endevco  7270A-2k  shock  accelerometers.  The 
excitation  was  provided  via  an  ICP  modal  sledgehammer.  The  frequency  response  functions  (FRF’s)  are 
measured  on  a  grid  using  a  roving  sensor,  fixed  input  approach.  Representative  data  is  shown  in  Figure  10(a.) 
and  (b.).  The  peaks  in  the  FRF  correspond  with  bending  modes  predicted  in  Figure  2. 
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Figure  10.  Schematic  and  representative  data  for  the  frequency  response  function  data  obtained  from  a  hanging 

plate  experiment. 

A  quick  estimate  of  the  plate  parameters  (i.e.,  dynamic  characteristics)  was  desired  from  the  experimental  data, 
notably  the  mode  shapes,  natural  frequencies  (ou),  and  structural  dampening  (£).  The  FRF’s  were  estimated 
using  a  simple  rational  fractional  polynomial  (RFP)  method  [14].  The  results  for  many  modes  simultaneously 
estimated  over  a  wide  frequency  range  are  shown  in  Figure  11  below.  Several  peaks  have  good  estimates,  yet 
several  others  are  not  accurately  captured  by  the  modal  estimates.  The  preliminary  modal  analysis  results 
presented  here  are  insufficient  as  a  characterization  of  the  mode  structure,  but  the  initial  indications  are  that  some 
of  the  bending  modes  are  being  predicted  accurately.  A  more  extensive  analysis  focused  on  the  hanging  plate 
using  rigorous  modal  analysis  techniques  is  ongoing  and  will  be  reported  in  a  future  paper. 
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Figure  11.  Estimates  on  the  mode  structure  of  the  hanging  plate  (inset  picture)  using  a  rational  fractional 

polynomial  parameter  estimation  algorithm. 


FUTURE  WORK 

As  noted  earlier,  several  observations  were  made  that  require  further  experimental  investigation.  New  data 
acquisition  hardware  is  being  used  to  capture  the  entire  impact  and  subsequent  decay  at  maximum  bandwidth. 
Additional  sensor  sets  are  also  being  examined;  for  example,  the  experimental  apparatus  is  being  modified  to 
include  integrated  quartz  crystals  to  create  effective  load  cells  (see  refs.  [15-17]).  An  ultra-high  bandwidth 
(8  GHz)  photonic  doppler  velocimeter  (PDV)  [18]  is  also  being  tested  for  surface-normal  velocity  measurements; 
this  provides  a  direct  indication  of  the  velocity  for  both  modal  estimation  and  computational  validation.  The  results 
from  these  studies  on  simple  structures  will  inform  the  development  and  implementation  of  higher  accuracy 
methods  for  numerical  dispersion  and  continue  with  studies  into  the  shock-response  behavior  of  complex 
structures  with  interfaces. 

SUMMARY 

The  dynamics  response  of  simple  bar  and  plate  structures  were  predicted  using  a  combination  of  analytic 
calculations  and  finite  element  simulations.  Experiments  to  reproduce  the  free-free  boundary  conditions  on  a  bar 
and  plate  were  described.  The  structures  were  excited  by  impact  from  impulse  hammers;  the  outputs  of  high 
bandwidth  strain  gages  and  accelerometers  were  monitored  to  determine  the  elastic  wave  propagation  in  the 
structures  as  well  as  the  corresponding  frequency-domain  characteristics:  power  spectral  density,  transfer 
functions,  and  coherence.  Leakage  was  found  to  significantly  reduce  sensor  coherence  over  a  broad  spectral 


band,  but  the  use  of  an  exponential  window  restored  the  coherence  in  the  original  signal.  Using  this  approach,  a 
broadband  coherence  was  obtained  up  to  approximately  10  kHz.  However,  resultant  broadening  of  the  modal 
peaks  alters  the  apparent  damping,  illustrating  the  difficulty  in  maintaining  coherence  in  a  laboratory  setting.  The 
initial  results  from  an  ongoing  effort  to  characterize  the  modal  response  of  a  hanging  plate  were  also  presented 
along  with  an  estimate  of  the  modal  parameters. 
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Research  Objectives 


Idesman,  Comput  Mech  (2007)  40:261-279 


Research  Challenge: 

•  Current  models  and 
numerical  methods  do 
not  agree  with 
experiments  at  high 
frequencies 

•  ^>10  kHz 

•  Can't  predict 
resonant  excitation 
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Linear,  Elastodynamic  Simulation  of  Bar  Impact 

Research  Objectives: 


•  Quantify  the  ability  of  "hydrocodes"  and  other  M&S  capabilities  to 
capture  full  spectrum  using  analytic  cases 


•  Validate  predictions  using  controlled  experiments 
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•  Perform  "simple"  validation 
experiments 

-  Simplified  "Hopkinson  bar" 

-  Resonating  Plate 

•  Apply  impulsive  loads 

-  Strain  rates  ~102-104s'1 

-  Force  »  100  N 

•  Use  high  bandwidth 
measurement  techniques 

-  f  >  10  MHz 
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Simplified  “Hopkinson  bar” 


Piezo¬ 

electric 

Stress 

Sensors 


Ti  bar 


Laser 
ppler 
velocimeter 


Accelero- 
Strain  gag§ieters 


Resonating  Plate 
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Work  by  Idesman 

-  Axisymmetric  FEA 

-  Measured  dispersion 

-  Compared  with  1-D  theory  (no  dispersion) 

Modal  analysis  using  LS-DYNA 


D=2R 

=2 


Free- free  bar  modes 
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Bar  I  mpact  Experiment 


Simplified  Bar  Experiment 

•  Short  length  to  accommodate 
high  mesh  densities 

•  Known  material  properties 

-  Ti  alloy  (6%  Al,  4%  V)  grade  5 

-  Length:  6  in. 

-  Diameter:  1  in. 

-  Weight:  0.754  lbs. 


Sensors 

•  Semiconductor  Strain  Gage 

-  5  mounted  with  epoxy 

-  Resistance  of  500  ft 

-  Gage  factor  of  140  ±  5%  at  75°F 

•  Piezoresistive  accelerometer 

-  Full  scale  of  2  kg 

-  Hard-mounted  on  end  of  bar 


Schematic  of  Strain  Gage  Attachment 


External  leads 
(enameled  wire) 


Solder  pad 
Gage  leads 
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Bar  I  mpact  Experiment 


•  Elastic  stress  wave 
introduced  with  impulse 
hammer 

-  Assume  center  point  impact 

-  Metal  tip  used 


Plastic  Rings 


Rubber 

Bands 
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Raised  30- 


Center  Point  Strain 

I  mpact  Gages  Accelero' 

meter 
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•  Analyzed  data  using  signal  processing  tools 

-  PSD 

-  Transfer  function 

-  Coherence 


•  Coherence  is  measure  of  correlation  between 
response  and  input  (stimulus)  energy 

-  Noise 

-  Causality  ^  |P,(/f 

 /•(/  )/'  (/  > 
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Transfer  Function  Estimate 


FtaquncpQQO 

Transfer  function  estimate  between  the  strain  gages  (data  is  filtered  at  50  kHz). 
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FtvqumcyQfiO 

Mean-squared  coherence  (data  is  filtered  at  50  kHz). 
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•  Long  sampling  time  required  for  low  frequency 

-  Very  long  record  lengths 

-  Competes  with  bandwidth 

•  Small  size  of  bar  creates  challenges  for  sensor  density 

•  M&S  requires  very  fine  meshes 

-  Results  in  very  long  simulation  times 

•  Experiment  data  much  higher  fidelity  and  coherence 


Even  simple  cases  are  computationally 
expensive  to  run  and  do  not  agree  with 
validation  experiments  under  impact  loads 
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Plate  I  mpact  M&S 


•  LS-DYNA  Simulations 

-  Modal  analysis 

-  Plate  impact 
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Free-free  plate  modes 
f1  =  22  Hz  f2  = 


27  Hz 
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Plate  I  mpact  Experiment 


•  Impact  Testing 

-  Free-Free 
Boundary  Conditions 

-  Eigenmodes 

-  2'  x  4'  Aluminum  Plate 


Sensors 

•  Force  hammer 


Accelerometers 

-  Endevco  shock  (2  k g) 

-  PCB  triaxial  modal  (500  g) 


•  PCB  piezoelectric  strain 
gages 
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Plate  I  mpact  Experiment 


X  Input 
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Data  Analysis 


•  Transfer  functions 

-  Force  hammer  to  strain 
gage  output 

•  Excellent  agreement  with 
modal  predictions 
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Transfer  Function  of  Strain  Gage  SN  3328 


Transfer  Function  of  Strain  Gage  Sn  3327 
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Data  Analysis 


•  Initial  modal 
parameter 
estimation 
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Parameter  Estimation  Using  Rational  Fractional  Polynomials 
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Results: 

•  Model-based  estimation  of 
dynamic  characteristics  from 
experimental  data 

•  Mode  shapes,  frequencies 
extracted  from  simulations 


Assessments: 

•  Modal  structures  agree 

-  Theoretical  (analytic) 

-  Simulations 

-  Experimental 
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Phase  and  Amplitude  Spectrum 


•  Bandwidth  vs.  record  length  (again) 

•  Mode  structure  is  observable  in  controlled  (lab)  conditions 

•  Initial  result:  models  accurately  predict  linear  modes 


Future  Work: 

•  Experiment  data  much  higher  fidelity  and  coherence 

•  Add  model-based  design  of  experiments  capability  for 
resonating  structures 

-  Impact-  or  explosively-driven 

-  Beams  or  plates 

•  Increase  loading  rate  and  amplitude 

-  Move  to  nonlinear  regime 
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The  dynamic  responses  of  simple  mechanical 
structures— bars  and  plates— were  modeled  and 
experimentally  measured 

At  low  loading  rates  (linear  problem)  mode 


structure  is  accurately  simulated 

•  Simulations  of  higher  loading  rates  are  indicating 
lower  effective  band  widths 

-  Evidence  in  coherence  of  experiment  vs.  M&S 

•  Future  work  will  focus  on  further  quantifying 
coherence  in  structures  and  developing  techniques 
to  increase  the  effective  bandwidth 
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